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1.  Introduction 


A thin  spherical  layer  of  an  incompressible,  inviscid  fluid 
which  is  held  on  the  surface  S of  a rotating  ball  by  gravitation 
can  be  taken  for  some  purposes  as  an  approximation  to  the  Earth's 
atmosphere.  An  investigation  of  the  two  dimensional  vortical 
motion  in  such  a layer  should  be  useful  for  the  understanding  of 
certain  observed  meteorological  phenomena.  For  example,  it  appears 
that  the  analysis  of  large  scale  closed  isobaric  systems  can  be 
based  on  a knowledge  of  the  paths  of  concentrated  vortices. 

This  report  presents  a study  of  the  vortical  motion  which  is 
due  to  the  existence  of  concentrated  vortices  (normal  to  S)  which 
are  confined  to  a polar  cap  and  subject  to  a boundary  condition 
along  a circle  of  latitude. 

We  assume  that  the  departure  of  the  free  outer  surface  of  the 
layer  from  an  equilibrium  position  is  small;  and  that  the  tangen- 
tial acceleration  is  negligible  compared  with  the  Coriolis  force. 

We  also  assume  that  the  variation  of  the  Coriolis  force  with  lati- 
tude can  be  neglected  In  other  words,  we  study  geostrophic  vor- 
tices on  a sphere  as  contrasted  wici  geostrophic  vortices  on  a 
plane.  The  latter  have  been  discussed  by  several  authors  in 
connection  with  the  motion  in  a rotating  planar  tangential  layer 
as  an  approximation  to  the  motion  in  a thin  layer  of  fluid  covering 
a rotating  ball.  References  can  be  found  in  the  paper  by 
Morikawa,  [ 1]  . 

The  main  result  is  the  system  of  nonlinear  equations  for  the 
paths  of  the  vortices.  These  are  used  to  study  the  linear  stabil- 
ity of  the  motion  of  n vortices  which  are  symmetrically  arranged 
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along  a circle  of  latitude  in  a cap  with  or  without  a polar  vortex. 
One  of  the  boundary  conditions  imposed  requires  the  velocity  to 
the  north  to  be  zero  along  the  circle  which  bounds  the  cap;  and 
the  other  requires  the  velocity  to  the  east  to  be  constant  along 
the  boundary.  The  investigation  is  restricted  to  analytical 
results  which  include  a formula  for  linear  stability  for  the  case 
in  which  the  polar  vortex  is  either  held  fixed  or  else  has  zero 
strength.  Numerical  computations,  which  are  necessary  for 
decisions  about  the  linear  and  nonlinear  stability  of  given  con- 
figurations of  vortices,  will  be  discussed  in  another  report.  See 
Leiva,  [2]. 
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2.  Equations  of  Motion 


Let  p denote  the  distance  of  a point  from  the  center  of  a 
ball  E of  large  radius  a which  rotates  with  constant  angular 
velocity  to  about  a polar  axis.  Let  and  0 denote  respectively 
the  longitude  and  the  colatitude  of  a point  on  the  rotating 
spherical  surface  S of  E.  Let  p = a and  p = a + h(<j>,0,t)  represent 
two  surfaces  which  contain  an  incompressible,  inviscid  fluid  which 
is  gravitationally  attracted  by  E.  Suppose  h is  small  compared 
with  a. 

The  velocity  of  a fluid  particle  relative  to  S is  defined  by 
the  components 


u = (p  sin  0)  = 

d0 

v = “P  “ 


w = 


tangential  component  toward  the  east, 
tangential  component  toward  the  north, 
radial  component. 


If  the  only  body  force  acting  is  that  due  to  the  gravitational 
potential  G of  E,  then  the  basic  hydrodynamicai  equations  which 
define  the  motion  of  the  fluid  relative  to  S are  the  continuity 


equation 


dw 

P 3p 


+ 2w  + 


1 

sin  0 


du  d(v  sin  0) 


0 


and  the  momentum  equations 


du  _ uv  cot  6 + uw  + 2uw  sin  0 - 2cjov  cos  0 
at  p p 


5oP 


sm 


5px 

d<{) 


5 


dv  + vw  + u cot_e  + 2cou  cos  e 
at  p p 

_ — - — _ 2cou  sin  b 
at  p p 


In  these  equations  the  differential  operator  with  respect  to  the 
time  means 

d = ^ , u d _ v 5 + a 
HT  "57  p sin  *3  p ^ ^p 

The  symbol  6q  denotes  the  constant  density  of  the  fluid;  and  p 
stands  for  the  pressure  in 

p p p 

6 p co  sin  6 
Pq  = P +6oG  “ “ 5 

which  can  be  referred  to  as  a modified  pressure. 

Since  h is  supposed  to  be  small  compared  with  the  large 
radius  a of  E;  and  since 

w(<l,  6,  a,  t ) = 0 


let  us  neglect  the  radial  velocity  and  the  radial  variation  of  u 
and  v.  Let  us  assume  that  G = gp  where  g is  a constant;  and  that 

the  centrifugal  effects  manifested  by  the  partial  derivatives  of 

, 2 2 . 2fl 
6 p a)  sin  b 

° can  be  ignored.  Let  us  also  assume  that  the  motion 

is  such  that  the  nonlinear  terms  in  the  tangential  momentum  equa- 
tions can  be  neglected;  and  that  the  radial  momentum  equation  can 
be  replaced  with  the  hydrostatic  law 


p(<M,p,  t)  = g^0(h  + a - p) 

which  satisfies  the  condition  that  the  pressure  is  zero  at  the  free 
surface  p = a+h. 
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Under  the  above  assumption  and  with  the  notation 


*u(  4>,  e,  t ) = u(cj>,e,a,  v)  , 


v(4>,e,  t)  = v($,  0,a,  t)  ; 


an  approximation  to  the  motion  is  determined  by  the  equations 


(2.1) 


dw  

"5p  a s 


1 r< 

TnT 


b(v  sin  0) 


= 0 ; 


(2.2) 


- 2o)v  cos  0 = 


g dh  . 
a sin  0 ’ 


(20) 


dv  . _ ~ g dh 

¥t  a«u  cos  9 = |¥ 


An  integration  of  (2.1)  gives 


w(M,a+h,t)  * a's'lnT 


d(v  sin  0) 


~W 


du 


The  kinematic  condition  at  the  free  surface  is 


dh 


w ( ij),  0,  a+h,  t ) - 


and  hence 
(2.4) 


dh 

cTE 


h 

a sin  0 


d(v  sin  0) 


~SS" 


du 

"5^ 


We  suppose  new  that  h(<{>,0,t)  is  always  close  to  the  constant  value 


h . Then  if  we  introduce 

o 


h - h 

rj ( 4>,  0,  t ) = — ft — 

o 


a linearization  of  (2.4)  yields 


\ 


3 . Geostrophlc  Vortices 


If  we  neglect  the  variation  of  the  Coriolis  force  with  lati- 
tude  and  take 

OJ  COS  6 = (D  COS  6^  * > 


then  the  last  equations  of  Section  2 reduce  to 

(5.1) 


^t  a sin  0 


d(v  sin  6)  _ du 

"5$ 


(5-2) 


(5.5) 


& - 


It  + Stv 


-w 


6h0  . 
a sin  5 * 


. gho  an 

a "5s 


Hereafter  we  confine  ourselves  to  a study  of  these  equations.  As 
win  be  explained  in  the  sequel,  they  lead  to  what  are  called  geo- 

strophic  vortices. 

The  elimination  of  tj  from  (5.2)  and  (5.5)  gives 


(5.4) 


The  quantity 


a 

■51 


a sin 


d(u  sin  6)  , bv 

58"  "5? 


= 2aJ1rit 


^ a sin  0 


5(u  sin  6 ) 5v 

^ 58"  "5$ 


is  the  radial  c 
we  have 

(3  *c  ) a sin  "0” 


omponent  of  vorticity;  and  by  integration  of  (5*4) 


5(u  sin  6)  5v 

” 58"  "5$ 


= 20)^  + C(<M>°)  ' 2a)iTi(<f»,e,°) 
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This  implies  that  we  can  use  (3*5)>  (3*2)  and  (3*3)  as  a basic  set 
of  equations  instead  of  (3*1);  (3*2)  and  (3»3). 

Equations  (3.2)  and  (3«3)  imply 


(3.6) 

v2  - 

d u , 

4<o^u 

_^o 

2all  ^ 

1 . 

i 

it? 

a 

) 

sin  6 

ata<p 

and 

(3.7) 

S2v  + 

4o)^v 

_^o 

' 2a)1 

an 

+ 

- 

i? 

a 

sin  8 

y 

J 

• 

If  these  are  used  to  eliminate  u and  v from  (3.5),  we  find 


(3.8) 


£ sln  6 & + ^ ^ ' Sh; 


^ + 


20 

a 2o>, 

U(<M,o)  -2co1TjU>,e,o)] 


a22co. 


The  theory  of  Laplace  transforms  can  be  used  to  show  that  the 
steady  state  solution  of  (3*8)  is  such  that 


(3.9) 


Voo1^'6’^  = 


where  ip  must  satisfy 


^ ~ sH7 


(3.10) 


sin  & 


sin  0 


sin  6 


aV(<M) 


The  function  ip  does  not  depend  on  t and  the  associated  time  inde- 
pendent velocity  components  are 


8 


(5.H) 


U = Lt->Co^,e,t)  = I H ; 


(5.12) 


v = Lt-oov(^^t} 


a sin  6 id 


Note  that  (5.10)  is  a consequence  of  the  equations  (3. 11),  (5*12) 
and  the  vorticity  equation 


(5.15) 


1 k(U  sin  6)  dV 
a sin  5 <5^ 


= 2a)l  Lt->ooT1+C*(^0) 


+ C*(4>,  e)  • 


In  other  words,  (3*10)  is  implied  by  (3*5)»‘  (5*2)  and  (3*3)  when 
we  ignore  inertial  forces. 

We  are  irterested  in  the  solution  of  (3*10)  when  the 
vorticity  £*(<j>,  0)  Is  constant  in  a small  neighborhood  of  a point 
((J)*,0*)  and  zero  elsewhere.  This  solution,  as  is  well  known,  can 
be  derived  from  the  solution  of  the  idealized  equation  which  comes 
from  (3.10)  when  £*(<[>, 0)  is  assumed  to  represent  a vortex  of 
strength  q concentrated  at  the  arbitrary  point  (<j>,0).  For  this 
reason  we  are  going  to  begin  with  the  mathematical  interpretation 
defined  by 


a2M.6((|)-<|)1)6(e-01) 

a2[C(M,0)  -2o>,T)($,e,0)  = a £ (<j>,0)  = 2 “ 

1 a sin  0, 


where  6 symbolizes  the  Dirac  delta  function.  Hence  the  basic  equa- 
tion to  be  solved  is 

.2  2 


(3.1*0 


sin  6 "3$ 


sin  0 


1 ^2  4a2o>2^  _ M.6(<t>-({>1)6(e-e1) 


sin20  d<|> 


sin  5- 


9 


I 


» 

> 

i 


* 


► > 


K 


i 

i 


\ 


L 


The  motion  in  a thin  planar  layer  of  fluid  tangential  to  the 
surface  of  the  Earth  is  often  used  as  an  approximation  to  the 
actual  motion  of  the  Earth's  atmosphere  in  the  neighborhood  of  the 


point  of 

tangency.  For  su 

rh  an  approximation  the  analogue 

(3.H); 

(3.12)  and  (3.1*0 

are 

(3.15) 

f 

s ; 

(3.16) 

< 

11 

dx  . 

and 

1 fdrU  , dV 

p.5  ( (J5- 4*!  )8  ( r_r!  ) 

(3.17) 

r {ST  + 

+ ri 

In  terms  of  the  polar  coordinates  (r,^),  these  show  that  x Inust 
satisfy 


(3.18) 


I a r 4*  + 
r 3r  75  ^ 


S;x 


He>(<t>-<t>1)6(r-r1) 


The  only  physically  admissible  solution  of  (3*18)  Is 


(3.19) 


M-  if 

X = - # Ko 


- — - J r2  + r^  - 2rr1  cos  (i-^) 


'gh 


where  K f ] denotes  the  zeroth  order  modified  Bessel  function  of 
o 

the  second  kind.  This  defines  what  is  called  a geostrophic  vortex 
The  motion  of  various  configurations  of  such  vortices  has  been 
studied  by  several  authors,  notably  G.  K.  Morikawa  [1],  [3]  whose 
papers  contain  detailed  explanations  and  references.  In  keeping 
with  what  seems  to  be  accepted  terminology  we  can  say  that  (3-1*0 


/ 


1 


/ 
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defines  a geostrophic  vortex  on  a sphere. 

Our  object  now  is  to  investigate  the  nature  of  the  solution 
of  (3. 14)  when  co  is  not  zero.  This  equation  can  be  written  as 


(3.20) 


sin 


a 

~u 


din  9 


— 4 v (v+1  )ip 
sin  9 


M-b  (<}>-<|>1)6(0-01) 
sin  6-^ 


where 

1,  2 2 

%).a 

(3.21)  v (v+1 ) = - -gp—  • 


The  solution  of  (3.20)  can  be  inferred  from  the  solution  of 


1 a 

sin  5 ~5& 


sin  6 ^ + v (v+l)^  = 


0 . 


This  equation  defines  the  Legendre  functions  of  degree  v.  Its 
general  solution  can  be  expressed  in  the  form 

t{/  = c1Pv(-cos  6)  + c2Pv(cos  6) 


when  v is  neither  zero  nor  an  integer  as  in  the  case  when  v is 
defined  by  (3.21)  from  which 


where 


v = - 


+ iq 


i = /-T  ; 


q.  = 


2 2 

1 I l6o^a 

&hn 


- 1 


and  / 0.  In  the  neighborhood  of  the  north  pole  Pv(cos  9)  is 
continuous  and  Pv(l)  = 1;  but  in  this  neighborhood 
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P (-cos  0)  = P , (-COS  0) 

V J-  I 4 ^ 


+iq 


= — cosh  qrr 

v 


CO 

I 


cos  qx  dx 


/2(cosh  x - cos  0) 


behaves  like 


, a \ 2 sin  vtt 

P (-cos  0)  ~ — • £n  S . 

v t r 


This  suggests  that  if  a vortex  is  concentrated  at  (^,0-^);  if  d-^ 
is  the  geodesic  distance  from  this  point  to  an  arbitrary  point 
( e ) on  S;  and  if  an  is  the  angle  o1  = dx/a,  then  the  fundamental 

solution  of  ("*  • ) is 


4 sin  vtt  Pv^"cos  al^  4 sin  mv  Pv 


-cos  0 cos  0- 


-sin  0 sin  6±  cos  ( )J 


It  can  now  be  verified  by  direct  substitution  that 


(3.22)  » - 4 sln~  Fv(-cos  °1> 

does  indeed  satisfy  (3.20).  Furthermore,  computations  with  the 
velocity  components 

U = 7t -A  P(-cos  o,  ) i 

4a  sin  vtt  o0  v 1 

V = 7'  ■•j  i I P ( — cos  Oi  ) } 

v 4a  sin  0 sin  vm  ^ v 1 

show  that  in  the  neighborhood  of  (^j  01),  (5*22)  defines  a vortical 
motion;  but  the  vortex  point  itself  remains  at  rest  - it  possesses 
no  autonomous  motion.  In  fact  it  can  be  shown  that 
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L 77 r= P (-cos  a-,  ) 

a*oo  4 sin  vt  v 1 

2oj 


u 


W Ko 


1 f_2  . 2 


r^  + r^  - 2rr^  cos 


Pence  we  can  say  that  (3.22)  represents  a geostrophic  vortex  on  a 
rotating  sphere  in  the  same  way  as  we  say  that 


X = 


/ 


? 2 

^ - 2rr1  cos 


(Mp 


represents  a geostrophic  vortex  on  a rotating  plane. 

For  an  arbitrary  distribution  of  n+1  vortices  on  the  sphere 
the  function  a stream  function  is 


n 


if  = 


sm  vir 


Li  ,P  (-cos 
vv 


where  a.  = d./a  and  d. 

J J J 

point  of  concentration 
point  (<{>,0).  That  is. 


is  the  geodesic  distance  on  S from  the 
of  the  j-th  vortex,  (ij,0j),  to  an  arbitrary 


cos  = cos  6 cos  0j  +sin  0 sin  0j  cos  ( 4>-4> j ) 


The  associated  velocity  field  is  given  by 


slnv,^{^jV-e°5  “Pi  ; 


V a sin  6 4a  sin  0 sin  vir  l u j 


{c^pv<-c0£ 


The  path  of  the  k-th  vortex  is  along  a stream  line, 
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(^(x,y)  = const),  and  its  velocity  is  equal  to  the  velocity  at 
^ue  oP  °^^er  vortices.  Therefore  since 

u = a sin  6 ; v = -a  ; 


the  equations  for  the  motion  of  (^  , 6 ) are 

K K 


d+k 


— 2 

4a  sin  6 sin  vtt 


w{^^Pv(-cos  °j)} 


JA 


e=e, 


and 


dek 

cTE 


— n 

4a  sin  6 sin  vtt 


°j)} 


J 

JA 


e=e,. 


If  the  motion  is  confined  to  a part  D of  the  sphere  S, 
boundary  conditions  must  be  imposed;  ard  it  is  not  then  possible, 
in  general,  to  express  as  a finite  sum  of  terms  of  the  type 

d i 

4 s In  vtt  Pv("cos  aj)*  When  D is  a cap  bounded  by  a circle  of 
latitude  we  have  a case  which  is  important  for  applications  to 
certain  meteorological  phenomena.  The  assumption  that  the 
Coriolis  force  is  independent  of  latitude  decreases  in  severity  as 
the  latitude  of  the  boundary  circle  of  the  cap  increases. 
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A 


4 . Ge ostrophic  Vortices  in  a Cap 


Suppose  that  n+1  concentrated  geostrophic  vortices  exist  in 
the  cap  which  covers  the  north  pole  and  is  bounded  by  the  circle 
of  colatitude  0 = X.  Suppose  that  along  this  boundary  circle  the 
velocity  ' the  north  is  zero.  In  order  to  study  the  motion  of 
the  vortices  we  need  the  stream  function  which  satisfies 


(4.1) 


sin 


Jt  (sin  0 ) ||  (<M  ) + — V"  “t!  + v (v+1^ 
in  R de  db  sin 6 

n jjl.6  (4>-4>.  )5  (0-0  . 

_ y — _3 3 -L 

Ifo 


sin  0 . 
J 


for 


0 <0,0.  < X ; 0 1 ^ ‘i'l  1 > 

— J 


with  the  boundary  conditions 


(4.2) 


V(fc,X)  = | = o 


or 


(4.3) 


x ) = c . 


The  solution  of  (4.1)  subject  to  (4.3)  can  be  expressed  in 
the  form 


n 


(4.4)  »(M)=bPv(cos  e)  + 4sl----  ^4jPv 


-COS  0 COS  0 . 

J 


-sin  0 sin  0^  cos  ($-<{)..) 


n oo 


4 sinvTr 


5S  £ 


C (m,X)F™(  cos  0)pm(cos  0 . )cos  m(|)-(j)  .) 
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where  pj(x)  is  the  spherical  harmonic  of  degree  v and  order  m. 
This  function  satisfies 


d I1-*2)  cl  ^<x)  + 


3x 


(v+i)  - -4 


l-xfc 


F^(x)  = 0 
v 


and  is  such  that 


f^D  = 


1 , m = 0 


0 , m / 0 . 


^(x)  can  also  be  expressed  in  terms  of  P^(x).  We  have 


(4.5) 


f^(x)  = (-lf(l-x2)m/2 


dmP  (x) 


dx 


m 


The  expansion 


P [-cos  6 cos  6 . - sin  6 sin  9 cos 


oo 


where 


and 


e F^(cos  6.)E^(-cos  9)  cos  m(  <|>-<t>,) 
m v j v J 


0<6<6<v  } <j>7  <t>j  real 


eo  = 1 J 

/ i 2r  (v  - m+3  ' 
em  - " r iv  + m-i  IT  ' 

P°(x)  = p„  (x)  '* 


m > 0 , 


shows  that  the  boundary  condition  (4.5)  is  satisfied  if  we  take 
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b 


and 


c 

P (cos  A) 

V 


C (m,A)  = -e 


m 


P^(-cos  A) 

_v_ 

F^cos  A) 

v 


It  can  be  shown  that  the  solution  (4.4)  is  unique;  and  that  it  is 
analytic  in  the  cap  minus  the  vortex  points. 

If  instead  of  (4.3)  we  impose  the  boundary  condition 


(4.6)  V/g(1tJ»^)  — c , 

which  means  that  the  velocity  to  the  east  is  kept  constant  along 
the  boundary  circle  of  the  cap,  then  (4.4)  again  provides  the  solu- 
tion if  we  now  take 


b = - 


sin  AP^  ( cos  A ) 


and 


(-cos  A) 


C (m,  A ) = e 

v m I*1  (cos  A) 


where  Pm'(x)  denotes  the  derivative  of  F^(x). 

v V 

The  series  which  appear  in  (4.4)  can  be  summed  if  the  cap  is 
the  northern  hemisphere.  For  this  case  A = tt/2  and 

cv(m  ’ £)  = + era  • 

Then  since 

oo 

(4.7)  YU  em^(cos  e)^(cos  6 ) cos  m(  ) 

rrR5  v v j j 

= P [cos  6 cos  6 - sin  6 sin  6 . cos  ($>-$>  .)] 

v j J J 
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0 < e,6j  < TT  ; e < 7T  ; <(>,  4>j  real  » 


we  see  that 


(4.8)  v,(<’»e)  °bpv(c°s  9)  * 4 sin  vir  ^ *3*' 


-cos  0 cos  0 


.1 


-sin  0 sin  0 cos  ( <fc-<J>j  ) 


n 


* 4 sin  v,r 


cos  0 cos  0 


-sin  0 sin  0j  cos 


J 


j}J 


where  for  the  boundary  condition 

= c 

the  upper  sign  is  to  be  taken  with 

v.  C 

b - TT^r  ' 

V 

while  for 

( 4 t = c 

the  lower  sign  is  to  be  taken  with 


Notice  that  P [-cos  0 cos  0 . - sin  0 sin0.  cos  is  singular  at 

($>,0)  = ( 4s  j'  0 j ) and  Pv[cos0cos0j-sin0sin0J  cos(*-fj)]  is  singu- 
lar at  ( 4>,  0 ) = (J>  ,tt  - 0 . ) - It  follows  from  the  form  of  (4.8)  that 
that  solution  can  be  constructed  from  the  fundamental  solution 


1 

4 sin  vtt 


P 

V 


[-cos  0 cos 


0 . - sin  0 sin  0_  cos  ( <|>-4>  1 )] 
J o J 
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by  using  a process  of  reflection  across  A = tt/2.  A reflection 
technique  for  the  construction  of  the  solution  for  a cap  other 
than  a hemisphere  is  unknown  when  v / 0. 

If  v = 0,  that  is  if  a)  = 0,  then  (4.1)  reduces  to 


(4.9) 


1 

sin 


a 

w 


sin  6 


"SB 


1 

TT7 

sm  6 


= 


n n16(<t>-<t>1  )&(6-6j  ) 
sin  F j 


This  equation  possesses  the  fundamental  solution 

6 - 6 


1 

W 


~ U C7  “1 

£n  a2  j tan2  | + tan2  - 2 tan  | tan  cos  ( )J 


= rr-  Re  £n  a 
2rr 


e i<f)  ei  ^j 
tan  e - tan  e J 


Let  us  also  remark  that  if  we  employ  the  transformation 

z = (a  tan  -^e1^  = x + iy 

then  the  problem  of  solving  (4.9)  subject  to  0 < 6 < A and  a 
boundary  condition  can  be  transformed  into  the  problem  of  solving 

+ t±  = q.6(x-x  )6(y-y  ) 

^2  Sy2  j=o  J J J 

subject  to  0 < J x2  + y^  < R and  a corresponding  boundary  condition. 
See  Peters,  [4j . The  latter  equation,  as  is  well  known,  governs 
the  motion  of  rectilinear  vortices;  and  for  many  cases  it  can  be 
solved  by  using  reflection  techniques. 

The  velocity  field  defined  by  the  stream  function  (4.4)  is 

given  by 
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20 


(4.10)  4a  sin  V7r*sin  ©k*ik  = _ijb  sin  V7r  sin  6kPv^COS 


n 


0^0 

jA 


[i  .P 

0 V 


-cos  6^  cos  6 j 


-sin  sin  6^  cos 


sin  6.  cos  0 . 
K-  J 


-cos  6.  sin  6 cos  (t|>  -(f>  . ) 
K 0 K J 


- sin  ev  2Z  V..  2Z  Cv(n.,x)^'  (cos  ek)f^<cos  6 ) cos  1.(^4 j ) > 
K j=0  3 m=G  v 


(4.11)  -4a2  s.  VTr-sin  6k*pk  = 


n 


*JPv 


-cos  6 cos  3 . 
K J 


-sin  0R  sin  6j  cos  ( ) 


. sin  ek  sin  0 sin  ( 4>k- 4> j 


■r,,rc>rf[  cos  ek  ><(  cos  6j  )m  sinn)(<J>k-<t>j  ) • 

3^5  J S5 


Suppose  now  that  one  vortex  of  strength  \xQ  is  at  the  north 
pole  and  that  n (n  > 2)  others  each  of  strength  u,  are  symmetric- 
ally  situated  on  a circle  of  latitude  within  the  cap.  In  other 
words,  suppose  that 

UvV  3 <4’0)  • 


u,  = u 
^o  o 


while 


(bye  ) b ((j-i)  ^,y) 


A 


j j 


j = 1,2,  ...,n 


J 


21 


where 


0 < 7 < A • 


For  this  configuration  the  velocity  to  the  north  of  the  polar 
vortex  is  determined  by 

-4a2  sin  VTT-eo  = n P^(-cos  7)  sin  7. sin  (J>0-$j) 


n 


' p-  Le. -*0 

K 


00  ^(cos  e.)  1 \ 

cv(m’xi  aln-g;-  <(coS  V”  Sln  <V*o)  • 

K. 


Using  (4.5)  we  find 

-4a2  sin  vt r*®0  (cos  7)  sin  7 sin  (<t,0“lt)j) 


C (1,  A)P^(l)Pv(cos  7)  sin  7 sin 


From  this  we  see  that 


0 


because 


The  velocity  to  the  north  of  any  one  of  the  vortices  on  the 
circle  of  colatitude  7 is  given  by 
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2 • 

-4a  sin  VTr.sin  -y»0. 


= p.  Y~  P1  [-cos27  - sin27  cos  (k-j)  ^]sin27  sin  (k-j)  & 


O'A 


n qd 


- d 


C (m.kHF^cos  7 )]2  m sin  m (k-j  ) ^ • 


J=i  m= 


In  this  expression  each  of  the  sums  with  respect  to  j has  the  form 


n 


^[cos  (k-j)  ~]  sin  (k-j)m-^- 


2tt 


This  sum  is  equal  to  zero  and  so  we  find  that 

K - 0 • 

The  velocity  to  the  east  of  any  one  of  the  vortices  on  the 
circle  6 = 7 is  given  by 

4a  sin  vtt-^  = -^b  sin  vttP  (cos  7) 


+ u [P*  (-cos  7)  -C  (0,  A)P*(cos  7)] 


o'-  v 

+ (j.  ) P^'[-cos27  - sin27  cos  j ““l2  cos  7 sin2J 


n 00 


- U 


C (m,  A )Pm*  (cos  7'An(cos  7)  cos  jm  ^ 


0=1  m= 


p 

4a  sin  vtt«Q  . 
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3l=1 


We  have  now  shown  that  if  at  time  t = 0 one  vortex  of 
strength  p is  at  the  north  pole  and  the  others  of  strength  p 
occupy  the  positions 

s ( ( J-l)  > 0 = 1,2,. ..,n  , 

then  the  polar  vortex  stays  at  the  north  pole  while  the  others 
remain  on  the  circle  of  colatitude  6 = -y  and  each  moves  along  it 
with  constant  angular  velocity  O.  The  strengths  pQ  and  p can  be 
chosen  so  that  0-0.  On  the  other  hand,  if  pQ  and  p are  pre- 
scribed it  may  be  possible  (depending  on  the  boundary  condition) 
to  choose  b and  y so  that  the  vortices  are  stationary. 
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5.  Linear  Stability  Equations 


The  motion  of  the  vox uices  described  in  the  last  section  is 
given  by 


(5.1) 


(*0*e0)  - ^o>°1 

( <tj»  e j ) - ((j-l)  ^+nt,-y)  , 1 = 1,2, ...,n  . 


In  order  to  investigate  the  stability  of  this  motion  consider  a 
second  sec  of  vortices  with  one  member  of  strength  nQ  situated  at 

(5.2)  ^o,6o)  = (ao+°t'f3o) 

and  with  the  others  of  strength  p.  situated  at 

(5*3)  (<|>j,0j)  = (aj  + (i”1)  'Tf  +Ot,-y+p  j ) , j = 1,2,  ...,n  . 


Suppose  that  at  t = 0 the  quantities  a^,  j = 1, 2,  ...,n;  and  (3^, 
j = 0, 1, 2, ...,n,  are  small.  If  these  quantities  remain  small  as 
time  goes  on  then  the  second  set  of  vortices  will  always  be  close 
to  the  first  set  (5.1).  A condition  for  linear  stability  can  be 
found  if  we  can  linearize  the  equations  which  result  when  the 
coordinates  of  the  second  set  of  vortices  are  substituted  in  the 
motion  equations  (4.10)  and  (4.11).  The  linearization  cannot  be 
performed  directly  because  aQ  need  not  be  small  in  order  to  have 
the  second  set  of  vortices  near  the  first  set.  However,  it  follows 
from  our  supposition  that 
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X 

CQ 

II 

cos  a 

a . 

0 

0 

0 

> 

0 

^ j — 1, 2,  • • • , n 

y0 

o 

cn 

ll 

sin  a 

0 

^0  V 

are  small  at  t = 0;  and  linearization  with  respect  to  these 
quantities  is  sufficient  for  the  purpose  of  developing  a linear 
stability  criterion. 

After  somewhat  lengthy  computations  the  linear  stability 
equations  turn  out  to  be: 


n-1 


2 . 2tt 


(5.4)  qj^-oy^  = (A0-B0^  sin  j — )yQ 


n 


+ °°w  S 


a.  COS  (J-l)  ^ + c'(7)  y—  p,  sin  (j-1)  2.  , 


n 


n o 


2tt 


J 


n 


n-1 


(5.5)  qi(y0  +ox0)  - -(A0  - B0  7—  cos2j  f))x0 


+ o0(-v)  2=  aj sln  (j-1)  *-coM  £ pj cos  ‘J-1’  ¥ - 


(5.6)  -qi  sin  7.pk  = ^2  Co(7)[xo  sin  (k-1)  ^ - yo  cos  (k-1) 


£hr- 


+ ^Vj  +^Bxj(“k-aj)  ’ 

jA  jA 


(5.7)  qx  sin  7-ak  = - ^ C^(7)[xo  cos  (k-1)  ~ + Y0  sin  (k-1) 


2tt- 


n 


JA 


aj  + ^ Ck/j  + Cf?k 

J'A 


where 
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^ = 


4a£ 


sm  V7T 

U. 


A = 


4b  sin  vtt 


H 


- py(u+^  p^(i)[cv(o,m -cv(i,x)p;an 


+ n[Cy (0, X )Py ( l)Py ( cos  -y) 

+ 2Cv(2,A)P^(l)P^(cos  7)  -P^(-C0S  7)  cos  7]  , 

Bo  = sin27[P^(-^os  7)  +4Cv(2,A)P^(1)P^(cos  7)]  , 

C (7)  = P'  (-cos  7)  sin  7 - Cv(l,A)P^(l)P^(cos  7)  sin  7 , 


*kJ 


2rr 


P^  (-cos  7kj)  sin  7 cos  7 sin  (k-j)  — 


_ 2tt 

+ (—cos  7kj)2  sin^7  cos  7 sin  (k-j)  ^ sin  (k-j)  ~ 

+ sin  7 .y2"  C (m,  ^ )Prn(  cos  y)^'(cos  y)m  sin  [m(k-j)  , 


2tt 


Bkj  = P^(-cos  7kJ)  sin  7 cos  (k-j)  — 
+ P^(-cos  7kJ)  sin1 7 sin2 (k-j) 


00 


C (m,  A ) [ F^(  cos  7)]2  m2  cos  [m(k-j)  -^-1  , 


2ij. 

= -P1  (-cos  7 . ) • [ sin27  + cos^7  cos  (k-j)  — ] 


Ckj 


7kj 


+ P^ (-cos  7kj).sinc;27  • sinH(k-j)  ^ 


sin27  .7“  C (m, A ) [ F^*  ( cos  7)]2  cos  [m(k-j)  , 


m=C 


27 


c 


4b  sin  vtt  -nil  / \ .2 

P (cos  7)  sin  7 

|1  v 1 ' ' 


+ — sin2y[P" (-cos  7)  +Cv(0,A)P^(cos  7)] 


V~  {p^(-cos  7 j ) cos  J -^-  + Pv(‘cos  7 j ) sin22y  • sin^j 
’1 

n 00  11  ^ ^ 

> ) C (m, A )f"  (cos  7 )P^ ( cos  7)  sin  7 cos  mj 

7 * V V 1 X 

OD  I OO 

Y~  cu(m>x)tpI!  (cos  7)]  sin  7 , 


P O Ojr 

cos  7kj  = cos  7 + sin  7 • cos  (k-j)  — , 

2 2 2tt 

cos  7^  = cos  7 + sin  7 cos  j — . 

The  linear  system  (5»*0-(5*7)  has  constant  coefficients. 
Therefore  the  investigation  of  the  stability  of  the  solution  of 

the  solution  of  the  system  can  be  started  by  substituting 

s t st  st  st 

XQ  = aQe  , o^.  = a^e  , yQ  = bQe  , 0 = bRe  for  the  variables. 

This  substitution  leads  to  the  characteristic  determinant  which 

theoretically  could  be  studied  in  the  usual  way;  but  it  is  clear 

that  the  required  analysis  for  the  general  case  would  be  too 

complicated  to  allow  the  deduction  of  a practicable  formula  for 

linear  stability.  If  a free  polar  vortex  is  involved  with  two  or 

more  other  vortices  it  is  best  to  subject  the  characteristic 

matrix  to  numerical  analysis.  C.  Leiva,  [2],  is  doing  this  and  he 

is  analyzing  the  results  for  the  case  in  which  the  cap  coincides 

with  the  northern  hemisphere. 

When  the  polar  vortex  is  either  absent  or  kept  fixed,  the 

system  (5.4)- (5.7)  simplifies  so  much  that  the  development  of  a 
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linear  stability  formula  becomes  feasible.  If  the  polar  vortex  is 

absent,  then  u.  = 0.  If  the  polar  vortex  is  present  but  con- 
o 

strained  to  remain  at  the  pole,  then  xq  = yQ  = 0.  Hence  for  either 
of  these  cases  our  system  reduces  to 


(5.8) 


(5-9) 


where 


n 


n 


-92\  ~ \/j  + Bkj(ak'aj)  ’ 

jA 


jA 


n 


n 


q2ak  = 


£ ^ + £ 


J 

jA 


ck/j+cpk 


jA 


q2  = Q1  Sin  7 = 


kar  sin  vtt  • sin  -y 


For  reasons  stated  above,  let  us  devote  the  remaining  analysis  in 
this  paper  to  the  system  (5*8),  (5*9)  and  its  implications. 

Summations  show  that  (5*8),  (5.9)  possess  two  invariants 
given  by 


(5.10) 

and 

(5.H) 

where 

D 
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r 


From  (5.10)  and  (5«ll)  we  obtain 


n 


k= 


\ nq2^1 


and 


n 


aR  = (C+Djn^t  +n^2 
where  ^ and  £2  are  small  quantities.  If  we  set 


gk  \ + q2^1 


a k = ak  + (C+D)^1t  + 


we  see  that 
(5.12) 


n 


^k  = 0 ' 


(5.13) 


n 


“k - 0 


and  we  find  that  (3^  and  must  again  satisfy  the  same  equations 


as  Pk  and  satisfy,  namely  the  system  (5.8),  (5*9).  In  matrix 


form,  we  have 


’ CK 

( N 

\ 

( ~ ■I 

bk 

(5.14) 

■ + q? 

• 

Uj 

'■V 

= 0 . 


Now  if  we  take 


e, 


V „St 
bk6 


~ ~ st 

ak = V 


then  (5.14)  reduces  to 


1 


JO 


(5.15) 


, ' + v 

1/ 

a. 

K1  '■aiJ 

= o . 


Here  the  submatrices  & , i]  and  £ are  right  circulant  nxn 
matrices.  It  can  be  shown  that  the  characteristic  roots  are  given 

by 


(5.16) 


q2slU)  = -A(|)  + /-BU)E(l)  1 

f H — 1,  2,  . . . , n ; 


q2s2(£)  = -A(£)  - /-BU)E(l') 


where 


A(i)  - 


n 


dA 


^ F' 
v 


n-1 


= -i 


P^  (-cos  7 .)  sin  7 cos  7 sin  j n 


2tt 

n 

. 2.7 r . . 2tt 


. + P"  (-cos  7 . ) *2  sin  7 cos  7 sin  j — sin  j — 

\ v J 

+ 7~  C (m,  A ) Pm( cos  7 ) F^1 ' ( cos  7 ) sin  7 . m sinjm  ~ 
v.  m=J  v v v< 


2tt 


• 2tt 
* sin  — 


which  is  a pure  imaginary;  while 


n 


BU)  = l B 


. 1 - 
) L 


i!(j-k)^T  1 


JA 


f ( — cos  7 )sin  7 cos  j n 


n-1 


(l-cosji^)  ^ + P”  (-cos  7j  )sin‘t7  sincj  n 

+ IZ  Cv(m,A)[^(cos7)]2n'2  cos  ~ J 


2 

1 - 

4 


2t r 
n 

2.  2tt 


2tt 


E(£)  = 


n 


c + ^Z."  c’--'e 


£(3-k)f 


j=- 

3 A 


kj 


= 4b  sln  VT-  P"(cos7)  sin2-, 


n 


+ sin^7[P"  (-cos  7)  + C (0,  *)P"  (cos  7)] 


(i 


n 1 ^ ^ (-cos  7 j ) [ c o s j ^L-  (sin27+  cos^y  cos  j^-)cos 


. 2tt 


2tt- 


3 


+ P"  (-cos 7 , ) sin22y  sin^j  ^ •(  1 + COS  j£ 


^ -v 
n 00 


. 2 


2tt 


2tt 


are  real. 

Since  a,  , 3,  are  supposed  to  be  small  at  t = 0,  the  quanti- 
k K 

ties  a,  , P must  be  small  at  t = 0.  We  define  the  motion  given  by 
K K- 

(5.1)  to  be  linearly  stable  if  c^,  remain  small  as  time 
increases.  Then,  in  order  to  have  linear  stability  it  is  necessary 
that  the  real  parts  of  the  roots  (5»l6)  be  less  than  or  equal  to 
zero.  Since  A(£)  is  a pure  imaginary,  it  follows  that 


(5.17)  B( £ )E( £ ) > 0 

is  a necessary  condition  for  linear  stability.  In  other  words, 
each  root  of  the  set  (5.16)  must  be  a pure  imaginary  number.  With 
this  condition  satisfied  we  still  need  to  consider  the  possible 
appearance  of  multiple  roots  in  the  set  (5«l6).  For  example  if 
£ = n then  A(n)  = 0;  B(n)  = 0;  and  consequently  there  is  at  least 
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a double  root  present.  Presumably  this  means  instability  since 
the  multiple  zero  root  indicates  solutions  of  (5*1^)  which  have 

IT  — 1 

the  form  t ' . However,  if  the  zero  root  is  just  a double  root  it 

can  be  shown  that  solutions  a which  in  part  depend  linearly 
on  t are  precluded  by  the  invariants  (5*12)  and  (5.13 )• 

We  notice  of  course  that 

+ ( C+D )£^t 

has  a part  which  depends  linearly  on  t.  This  is  a reflection  of 
the  fact  that  if  the  vortices  (5*1)  are  slightly  displaced,  then 
the  angular  velocity  0(7)  will  change  by  a small  amount.  The 
quantity  (C+D)  is  equal  to 

C+D  = ka.  sin  vtt  • sin  7 . 

Double  roots  other  than  a double  zero  root  may  imply 
instability . 

If  u.  =0  the  expectation  is  that  at  least  one  domain  of 
o 

linear  stability  is  the  cap: 

0_^7<As<A. 

In  fact,  if  ur  has  the  same  sign  as  11,  and  the  polar  vortex  is 
kept  fixed,  then  an  investigation  of  the  signs  of 

B ( £ ) , E(£),  £ = 1,  • • • n-1 

shows  that  the  product  B(£)E(£)  is  posits  e (no  matter  what  n is) 
provided  7 is  sufficiently  small.  It  is  interesting  to  note  that 


I 
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for  the  case  of  a ring  of  Bessel  vortices  the  linear  stability 
condition  is  not  satisfied  for  the  neighborhood  of  the  origin  if 
n The  boundary  angle  0 = Ag  can  be  approximately  determined 

by  imputing  B(£)  and  E(£)  for  various  values  of  -y.  For  any  given 
case  this  would  require  considerable  numerical  analysis.  As  an 
example  of  what  is  involved  we  turn  to  the  case  of  three  vortices 
on  a circle  of  latitude. 
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6.  Three  Vortices  on  a Circle  of  Latitude  in  a Cap,  Polar  Vortex 


Fixed  or  Absent 

Our  intention  is  to  devote  the  rest  of  this  paper  to  a 
discussion  of  the  case  of  three  vortices  of  equal  strength  q 
symmetrically  situated  on  a circle  of  colatitude  0 <_  6 < A.  When 
it  is  admitted,  a polar  vortex  of  strength  p.  is  assumed  to  be 
held  fixed  at  the  north  pole.  We  are  interested  in  the  linear 
stability  of  motion  of  the  vortices  whose  positions  at  any  time  t 
are  given  by 

(♦o’V  s <*o’0)  * 

(♦0>8j)  = (O-l)  ^+0t,7)  , J = 1,2,3  . 


Our  primary  object  here  is  not  to  present  an  extended  numerical 
analysis,*  but  rather  to  record  the  quar  ities  which  need  to  be 
evaluated  for  such  an  analysis. 

The  angular  velocity  of  the  vortices  in  the  equilibrium 
position  is  given  by 

2 

(6.1)  la.  .sAn  V1r  -n  = djh-iin  " ■ P' (cos  7) 

+ [P^t-cos  7)  - Cv(0,A)pJ(cos  7)] 


i 2 12 

+ 3Pv(-cos  7 +?y  sin  7)»cos  7 


co 


m= 


Cv(m,A)P^  (cos  7 )F^(  cos  7MI  + 2 cos  m c-y)  . 


For  the  present  case  of  three  vortices  we  find 
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B(  1 ) = B( 2 ) = B , 
E(l)  = E( 2 ) = E 


(6.2)  B = - ^ Pv(-cos2-y  sin27)  sin*7 


+ | P^(-cos27+^  sin2-y ) sin\ 


C (m,  A)[P^(cos  7)]2  m2  cos  m 


(6.J  ) E = 4b-sbn  . p^(cos  7)  sin27 


+ -y  [Pj(-cos  7)  +Cv(0,A)P^(cos  7)]  sin  7 


3 p. 
2 Pv 


P^(-oos27+-^  sin27)  cos27 


+ p;(-cosS  +-J  sinc7 ) Sinc27 


2 , 1 _ . _2 


+ sin27  » ) (m,  A)F^(cos  7)P^  (cos  7 ) * ( 1 + 2 cos  m ^-) 

+ sin27  • ) C (m,  A ) [ P^1  (cos  7 )] 2 ( 1 - cos  m ^L)  . 


For  & = 3,  B(3)  = E(3)  = 0.  This  means  a double  zero  root  in  the 
set  (5.16)  for  three  vortices.  There  are  no  other  multiple  roots 
and  hence  the  linear  stability  criterion  is 


B(i)EU)  > 0 , i = 1,2 
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(6.5) 


BE  > 0 . 


If  (6.5)  holds  then  the  component  periods  of  the  motion  about 
the  equilibrium  position  are  given  by 


sUi  _ A(l)  1 1,/BE  i 2 

2rri  SFTq^ 

where 

(6.6)  iA(  1 ) =|  P^(-cos27  + | sin27)  sin  7 cos  7 

+ P"  ( -cos2^  t-j  sin2-y ) sin^^  cos  7 


and 


Cv(m,A)P^(cos  7)^ 


(cos  7 ) m sin  m 


2tt 

T 


A(2)  = -A( 1 ) . 

Recall  that  if  the  boundary  condition  is 

^(<t>»A)  ~ c 


then 


13  p (cos  a)  ’ 

P (-cos  A) 

CV(0,A)  = " P^(cos  xr  ’ 


Cv (m, A) 


-2(-l) 


m r (v  - m+1) 
r ( v + m+l) 


Pm(-cos  A) 

V 

F^(cos  A) 


m ^ 0 . 


If  the  boundary  condition  is 


(<!>»*)  = c 
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then 


b 


sinAP'(coF  A) 


Cv(m,  X) 


P1  ( -cos  X ) 

c (o,a)  = 4- > 

v Pv(cos  A) 

n(  1 ,m  r(v  - m+1)  . ^ ('cos 

( } ^ ^7’ 


m ^ 0 


If  the  cap  coincides  with  the  sphere  then 


b = 0 ; Cv(m,Tr)  = 0 . 


Some  ideas  about  the  effects  of  imposing  a boundary  condition  for 
an  arbitrary  cap  can  be  drawn  from  a comparison  of  the  motion  con- 
fined to  the  northern  hemisphere  with  the  free  motion  on  the  whole 
sphere.  If  the  cap  is  the  northern  hemisphere  then 

°v(m’|)  - Tem  , 


= 1 


c _ p/  1 r(v  -m-l) 
ein  ' r (v  + m+1 ) 


m > 0 


and  the  series  which  appear  in  the  above  expressions  can  be  summed 
by  using  (4.7).  We  find  that  for  A = ^ : 
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\ 


6.7) 


4a_sin_^r  -4b  sin  vtt  . p.(cos  ) 

u u U Vv  r' 

+ [p;(-cos  7)  ± p;(cos  7)] 

P 1 2 

+ 3P^(-cos  7+7;  sin  7)  cos  7 

1 2 12 

± [ 2P^  ( cos  27)  + Pv(cos  7 +-^  sin  7)]  cos  7 

(6.8)  iA(  1 ) = | P^(-cos27  +|  sin27)  sin  7 cos  7 

+ P^(-cos27  sin27 ) sin^7  cos  7 

2 1 2 

+ [P^  (cos  7 + ^ sin  7) 

- ^ P^  cos27+^  sin27)  sin27]  | sin  7 cos  7 

A(2)  = -A(l)  , 

(6.9)  B = - | P^(-cos27+i  sin27)  sin27 

+ | P^(-cos27+^  sin27)  sin\ 

- [ P^  ( cos27  + sin  7) 

- \ P"(cos27+-^  sin27)  sin27]  ^ sin27  , 


(6.10)  E = — ~~  P^(cos  7)  sin27 


+ ^ [P;(-cos  7)  *P^(cos  -y)]  sin27 


\ P^  (-cos^7 +•£  sin^7)  cos^7 


2 . 1 _ . 2 


+ (-cos^y  + 4j  sin^y)  sin^P-y 


2 . 1 .2  , _ . 2, 


f P^  (cos  2^)2  sin7  + P^  ( cos  27)2  sin^27 


- P^(cos27+^  sin27  ) ( 2 - cos27) 
+ P^  ( cos27  + sin27)  • -£■  sin227 


If  the  boundary  condition  is 


^(i  » 5)  = c 


the  upper  sign  is  to  be  taken  with 


b = 


If  the  boundary  condition  is 


V'gO t1  > - c 


the  lower  sign  is  to  be  taken  with 


b = - 


pToT 

v 


A numerical  analysis  for  the  caps 
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0 < 6 < v 


0 ^ 9 < 


will  be  presented  in  a separate  report  by  C.  Leiva  [2]. 


For  the  arbitrary  cap 

0 <_  6 < A < 

each  series  which  appears  In  the  expressions  (6.1),  (6.2)  and  (6.6) 

can  be  approximated  by  the  sum  of  its  first  n terms.  Various 
methods  can  be  used  to  estimate  the  error  involved.  One  method  is 

to  use  the  formula 

o -m/ 2 1 

v , , xin  r(v  -m+1)  . = t1"* \ / (t-x)m-  P.(t)dt 

(6.11)  (-1)  p(v  +m+rj  V ; r(m)  J v 

X 

Pv(t)  ,l-x,m/2 

” flrniTT  ’ TFx' 
x < I < l ; m^l1 


and  the  formula 


(6.12) 


,m  r (v  - ro+1 , 
' r (v  +m+i 


t-C'wi 


Pv(t)  ,1-X,m/P  1_ 

-rw  lr+v  i.x 


in  conjunction  with 

(6.15) 


t xm  r(v  -m+1)  = 

^ ' p (v  + m+1 ) 


m > 2 


cosh  qm 

irr^tojp”'  (0) 


For  example,  consider  the  series  which  appears  in  the  expression 


for  O,  namely 


oo  ml  _ ^ 

= - cv (cos  •Y)Prn(cos  y)(l  + 2 cos  m ) . 
m=0  v v 


Let  the  boundary  condition  be 


so  that 


Cv(m,A) 


^(<{bA)  = c 


?v(-cos  A) 
Py(cos  A) 


■ 2 ( - 1 )m 


r(v  - m+l)  ^(-cos  *) 

FT7~F5TTT  ( cos  X)  ’ 


m ^ 0 . 


If  Sg,  the  sum  of  the  first  two  terms,  is  used  to  approximate 
we  have 

, .V-  - 6 cosh  g.  *2.  M Pf'(c°s7)lf(cos7) 

2 T 5 Pm(cos  X)  Pim'(0)P}m(0) 


Then  from  the  formulas  (6.11)  and  (6.12)  we  find 


P (cos  A ) [ P (0)1 2 m= 


00  (cot2  » tan2  l)3m 


* b !J=J  < Pv<-cos  M[Pv(cos  7)]2  ]—  (cot2  ! tan2  5f) 


2 A .2  -vx^m 


Since 


we  obtain 


C = cot  ^ tan  < 1 
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.12 


Fv(cos  A)[Pv(0)r  1 -C 


^(s2  ~YZ  ) 

b cosh  qTr 


P (-cos  A ) [ P (cos  y)]^12 

V V 

1-^ 


This  shows  that  YZ  - 00  as  C - !»  which  means  that  ° iS  large 
when  the  vortices  on  the  latitude  circle  are  close  to  the  boundary 

circle.  It  also  shows  that 


S2" 


— 0 


as  7 - 0.  A similar  analysis  can  be  applied  to  the  infinite  sums 
in  (6.2),  (6.3)  and  (6.6). 
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